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ABSTRACT 



The relatively recent insight that energy input from supermassive black holes 
(BHs) can have a substantial effect on the star formation rates (SFRs) of galaxies 
motivates us to examine the effects of BH feedback on the scale of galaxy groups. At 
present, groups contain most of the galaxies and a significant fraction of the overall 
baryon content of the universe and, along with massive clusters, they represent the 
\ s 2^ ! only systems for which it is possible to measure both the stellar and gaseous bary- 

■ onic components directly. To explore the effects of BH feedback on groups, we analyse 

two high resolution cosmological hydrodynamic simulations from the Overwhelmingly 
Large Simulations (OWLS) project. While both include galactic winds driven by su- 
pernovae, only one of the models includes feedback from accreting BHs. We compare 
£^ | the properties of the simulated galaxy groups to a wide range of observational data, 

including the entropy and temperature profiles of the intragroup medium, hot gas 
mass fractions, the luminosity— temperature and mass— temperature scaling relations, 
the K-band luminosity of the group and its central brightest galaxy (CBG), star for- 
• mation rates and ages of the CBG, and gas- and stellar-phase metallicities. Both runs 

C$ ' yield entropy distributions similar to the data, while the run without AGN feedback 

yields highly peaked temperature profiles, in discord with the observations. Energy 
input from supermassive BHs significantly reduces the gas mass fractions of galaxy 
groups with masses less than a few times 10 14 M Q , yielding a gas mass fraction and 
X-ray luminosity scaling with system temperature that is in excellent agreement with 
the data, although the detailed scatter in the L — T relation is not quite correct. The 
run without AGN feedback suffers from the well known overcooling problem — the 
resulting stellar mass fractions are several times larger than observed and present-day 
cooling flows operate uninhibitedly. By contrast, the run that includes BH feedback 
yields stellar mass fractions, SFRs, and stellar age distributions in excellent agreement 
with current estimates, thus resolving the long standing 'cooling crisis' of simulations 
on the scale of groups. Both runs yield very similar gas-phase metal abundance profiles 
that match X-ray measurements, but they predict very different stellar metallicities. 
Based on the above, galaxy groups provide a compelling case that feedback from 
supermassive BHs is a crucial ingredient in the formation of massive galaxies. 

Key words: galaxies: formation — galaxies: clusters: general — cooling flows — 
cosmology: theory — X-rays: galaxies: clusters — intergalactic medium 
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1 INTRODUCTION 

Galaxy formation is one of the key unsolved problems in 
modern astrophysics. Traditionally, studies of galaxy for- 
mation have focused on the observable stellar properties of 
galaxies. However, a large fraction of the baryonic mass of 
massive galaxies is believed to be in a hot diffuse form and 
many of the processes that regulate the formation of stars 
do so by influencing the properties of this hot gas. Thus, a 
complete view of galaxy formation necessarily incorporates 
both the stars and the hot gas and an understanding of the 
processes by which these phases interact. 

At present, galaxy groups and clusters represent the 
only systems in the universe for which it is possible to mea- 
sure the properties of both the stellar and gaseous baryonic 
components out to a significant fraction of the halo viral ra- 
dius. This is possible because the gaseous component is com- 
pressed and heated to very high temperatures (T ~ 10 6-8 
K) where this material (the intragroup or intracluster me- 
dia) becomes visible at X-ray wavelengths. While clusters 
have higher X-ray surface brightnesses than groups, and the 
derived physical properties of the gas are therefore presum- 
ably more robust, clusters represent very rare and extreme 
environments. Since the efficiencies of many of the environ- 
mental processes that are thought to alter galaxy forma- 
tion/evolution (e.g., ram pressure stripping, strangulation, 
tidal stripping, galaxy harassment) increase rapidly with the 
mass of the main halo, one expects that galaxies in clusters 
have been subjected to far more violent interactions than 
those which reside in the lower density and more typical en- 
vironment of groups. It is therefore possible that a picture 
of galaxy formation derived in large part from observations 
of clusters is misleading. 

Galaxy groups, on the other hand, contain a large frac- 
tion of present day galaxies and a significant fraction of the 
overall universal baryon budget (e.g., Mulchaey 2000). They 
are therefore more representative hosts than clusters. In ad- 
dition, the effects of feedback from supernovae (SNe) and 
supermassive black holes (BHs) on the hot gas, which is 
thought to be crucial in shaping the properties of galaxies, 
are expected to be much more obvious (and therefore per- 
haps easier to quantify) on the scale of groups, where the 
energy input associated with these sources is comparable 
to the binding energies of these systems. And while what 
could be said about the detailed properties of the hot gas 
was severely limited in the past (with data from the previous 
generation of X-ray satellites), new high-quality data from 
the Chandra and XMM-Newton X-ray telescopes has yielded 
a much more detailed picture of the thermodynamic state 
of the gas in groups. For these reasons, galaxy groups prob- 
ably represent the best objects for studying gas and stars 
simultaneously. 

Much progress has also been made in the theoretical 
modeling of the formation and evolution of groups and clus- 
ters. For example, the past two decades have seen the de- 
velopment of sophisticated cosmological hydrodynamic sim- 
ulation codes capable of simulating the properties of the hot 
gas in groups and clusters ab initio with millions to hundreds 
of millions of particles/cells (e.g., Evrard 1990; Thomas & 
Couchman 1992; Navarro et al. 1995; Bryan & Norman 1998; 
Kravtsov et al. 2000; Borgani et al. 2001; Springel et al. 
2001; Kay et al. 2002). Broadly speaking, these simulations 



have had varying degrees of success in reproducing the ther- 
modynamic properties of the hot gas but, in general, those 
simulations that included the effects of radiative cooling ap- 
pear to have formed far too many stars relative to what is 
observed in local groups and clusters (see Balogh et al. 2001 
for further discussion of this 'cooling crisis'). This may be 
because, until very recently, such simulations did not include 
feedback from supermassive BHs. In recent years, a growing 
number of authors have argued, based on the results semi- 
analytic models of galaxy formation as well as hydrodynamic 
simulations (usually idealised, as opposed to cosmological), 
that feedback from supermassive BHs plays a crucial role in 
regulating the star formation rates of massive galaxies (e.g., 
Springel et al. 2005, hereafter SDH05; Bower et al. 2006, 
2008; Croton et al. 2006; Hopkins et al. 2006; Somerville et 
al. 2008) and suppressing the onset of catastrophic cooling 
of the hot gas in groups and clusters (e.g,. Churazov et al. 
2001; Brighenti & Mathews 2003; Dalla Vecchia et al. 2004; 
Ruszkowski et al. 2004; Omma et al. 2004). 

The insight that AGN feedback is important in the for- 
mation and evolution of massive systems motivated SDH05 
(see also Thacker et al. 2006) to develop a novel scheme 
for following the growth of supermassive BHs and feedback 
from AGN self-consistently in cosmological smoothed parti- 
cle hydrodynamics (SPH) simulations. Using the model of 
SDH05 as a backbone, Sijacki & Springel (2006), Sijacki 
et al. (2007), Puchwein et al. (2008), Bhattacharya et al. 
(2009), and Fabjan et al. (2009) examined the effects of AGN 
feedback on groups and clusters. These authors indeed find 
that including self-consistent AGN feedback yields reduced 
stellar mass fractions and star formation rates while repro- 
ducing a number of key hot gas observable properties as well. 
For example, Puchwein et al. (2008) found that including 
energy input from supermassive BHs enables their simula- 
tions to reproduce the luminosity-temperature and gas mass 
fraction-temperature relations of groups and clusters. Fab- 
jan et al. (2009) found similarly good matches to these re- 
lations, while demonstrating also that the simulations yield 
metallicity and temperature profiles in reasonable agreement 
with observations of galaxy groups (although the predicted 
stellar fractions still appear to be somewhat too high and 
the entropy-temperature scaling is shallower than observed). 

In the present study, we examine the effects of AGN 
feedback on galaxy groups using a novel set of simulations. 
In particular, we extract two runs from the Overwhelmingly 
Large Simulations (OWLS, for short; Schaye et al. 2009) and 
attempt to ascertain whether or not feedback from super- 
massive BHs is a necessary ingredient for explaining the 
properties of observed galaxy groups. The simulations we 
analyse have very different implementations for radiative 
cooling, star formation, chemo dynamics, and feedback from 
supernova and AGN from those of SDH05 and the studies 

1 Note that Sijacki et al. (2007), Puchwein et al. (2008), and Fab- 
jan et al. (2009) modified the original implementation of SDH05 
to include both 'quasar' and 'radio' feedback modes. The former 
is characterised by isotropic thermal heating of the gas when ac- 
cretion rates are high, while the latter, which is imposed when 
accretion rates are low, is characterised by thermal heating of 
bi-polar spherical regions offset from the galaxy hosting the BH, 
which arc meant to mimic 'bubbles' observed in many nearby 
clusters. 
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mentioned above which are based upon that model. This 
means a comparison of the findings can be used to check the 
robustness of the conclusions to differences in the param- 
eterisations of various subgrid physics as well as inclusion 
of new physics. In addition, because the OWLS simulations 
are performed in large periodic boxes (as opposed to being 
re-simulations using 'zoomed initial conditions') they yield 
much larger samples of groups for analysis than used previ- 
ously for simulations with AGN feedback, albeit at lower 
mass resolution (e.g., each OWLS 100 h _1 Mpc box run 
yields « 200 groups with masses greater than 10 13 M , 
whereas the studies based upon re-simulations typically only 
analyse a handful of systems). This is potentially impor- 
tant, as observed groups and clusters show a large degree 
of system-to-system scatter at fixed mass in the properties 
of the hot gas (e.g., Fabian et al. 1994; Markevitch 1998; 
McCarthy et al. 2004). We compare the OWLS simulations 
to a much wider range of X-ray and optical observations of 
galaxy groups than previously considered. This allows for 
a more thorough assessment of the successes/failures of the 
models. Finally, the large suite of OWLS runs also allows 
us to explore the impact of additional sub-grid physics, not 
just for impact of AGN feedback. For example, the effects 
of changing the supernova feedback prescription, switching 
off metal-line, and changing the stellar initial mass function. 
These variations will be explored in a forthcoming paper. 

The present paper is organised as follows. In Section 2 
we provide a brief description of the OWLS runs examined in 
this study and how we select and analyse the galaxy groups 
that form in these simulations. In Section 3 we make detailed 
like- with- like comparisons with a wide range of observational 
data of both the hot gas and stellar baryonic components. 
In Section 4 we summarise and discuss our findings. We 
adopt a flat ACDM cosmology with h = 0.73, Q b = 0.0418, 
fi m = 0.238, and cr 8 = 0.74, which are the maximum- 
likelihood values from fitting to the 3-year WMAP CMB 
temperature anisotropy data (Spergel et al. 2007). When 
comparing observed and simulated metal abundances, we 
consistently adopt the solar abundances of Asplund et al. 
(2005). 

Finally, we note this is the first in a series of papers 
investigating galaxy groups in the OWLS suite of simula- 
tions. In Paper II we examine in detail the thermal history of 
the intragroup medium and present a simple physical model 
for the excess entropy seen at intermediate/large radii in 
groups. In Paper III we explore the effect of varying the 
parameters of the key subgrid physics models (e.g., stellar 
IMF, chemodynamics, supernova and BH feedback) on the 
properties of galaxy groups. In Paper IV we investigate the 
build up of metals in the intragroup medium and the origin 
of the large scale gradients. 



2 GALAXY GROUPS IN OWLS 

An in-depth presentation of the OWLS project can be found 
in Schaye et al. (2009). For the present study, we shall there- 
fore give only a brief description of the simulations and our 
analysis methods. In Paper III in our series on galaxy groups, 
we will provide a detailed description of our analysis meth- 
ods, including how we find and select the simulated galaxy 



groups and then produce and analyse mock observations of 
them. 

2.1 Simulation characteristics 

The OWLS project has been undertaken with the goal of 
exploring galaxy formation and its sensitivity to both re- 
solvable and 'sub-grid' physics in fully self-consistent cosmo- 
logical hydrodynamic simulations. The approach is, starting 
from identical initial conditions, to vary the key parame- 
ters of the physics modules (e.g., such as the adopted stellar 
initial mass function (IMF), the rate at which black holes ac- 
crete mass, the wind velocity and mass-loading of supernova 
driven winds, switching on/off metal-dependent radiative 
cooling, switching on/off mass/energy transfer from SNIa 
and AGB stars, etc.) to see what effects they have on the re- 
sulting population of galaxies, or in this case galaxy groups. 
Comparisons between the various runs should allow one to 
isolate which physical processes are most important for, for 
example, establishing the shape of the stellar mass function 
or the entropy of the intragroup medium. 

In Papers II and III in our series on galaxy groups, we 
will present detailed comparisons between the various OWLS 
runs and focus on the physics of the hot gas in galaxy groups. 
In the present study, however, we select just a subset of two 
OWLS runs; the reference run (hereafter REF) and the de- 
fault AGN feedback run (hereafter AGN). These two runs 
are identical in every way, except that the latter incorpo- 
rates a prescription for BH growth and AGN feedback, as de- 
scribed in Booth & Schaye (2009) (hereafter BS09), whereas 
the former does not. The goal of the present study, there- 
fore, is to ascertain whether including feedback from AGN 
is required to explain the observable properties of galaxy 
groups. 

The OWLS runs are initialised from identical initial con- 
ditions; i.e., from a ACDM cosmological density field (with 
a power spectrum computed using CMBFAST) in a peri- 
odic box of 100ft -1 Mpc on a side at z — 127. In generat- 
ing this density field, the various relevant cosmological pa- 
rameters (0,1,, 0, m , Qa, h, as, and n s ) were set equal to 
their maximum-likelihood values found from the analysis of 
the 3-year WMAP CMB data (Spergel et al. 2007). The 
simulations were evolved to z = using the TreePM-SPH 
code GADGET-3 (last described in Springel 2005) and we 
extract from the final snapshot (z = 0) all galaxy groups 
for which M200 ^ 10 13 M (where M200 is the total mass 
within a radius that encloses a mean density of 200 times 
the present-day critical density of the universe). The simu- 
lations use 2 x 512 3 particles, yielding a mass resolution for 
gas and dark matter particles m gas « 8.65 x 10 7 /i _1 Mq and 
m dm « 4.06 x lO 8 ^ 1 M s , respectively. Thus, a ~ 10 14 /i _1 
M Q galaxy group is resolved with ~ 10 5 gas and dark matter 
particles. As we demonstrate in the Appendix, this resolu- 
tion is sufficient to robustly predict the properties of galaxy 
groups in our simulations. 

The GADGET-3 code has been substantially modified 
to incorporate new baryonic physics. Radiative cooling rates 
for the gas are computed on an element-by-element basis by 
interpolating within pre-computed tables (generated with 
CLOUDY) that contain cooling rates as a function of den- 
sity, temperature, and redshift calculated in the presence of 
the cosmic microwave background and photoionisation from 



4 I. G. McCarthy et al. 



a Haardt & Madau (2001) ionising UV/X-Ray background 
(see Wiersma et al. 2009a). Star formation is tracked in the 
simulations following the prescription of Schaye & Dalla Vec- 
chia (2008). Gas with densities exceeding the critical den- 
sity for the onset of the thermo-gravitational instability is 
expected to be multiphase and to form stars (Schaye 2004). 
Because the simulations lack both the physics and the reso- 
lution to model the cold interstellar gas phase, an effective 
equation of state (EOS) is imposed with pressure P tx p 4 ^ 3 
for densities % > n* where n, = 0.1 cm~ 3 . As described 
in Schaye & Dalla Vecchia (2008), gas on the effective EOS 
is allowed to form stars at a pressure-dependent rate that 
reproduces the observed Kennicutt-Schmidt law (Kennicutt 
1998) by construction. The timed release of individual ele- 
ments ("metals") by both massive (Type II SNe and stel- 
lar winds) and intermediate mass stars (Type la SNe and 
asymptotic giant branch stars) is included following the pre- 
scription of Wiersma et al. (2009b). A set of 11 individual 
elements are followed in these simulations (H, He, C, Ca, 
N, O, Ne, Mg, S, Si, Fe), which represent all the important 
species for computing radiative cooling rates. 

In the REF run, feedback from supernovae is incor- 
porated using the kinetic wind model of Dalla Vecchia & 
Schaye (2008) with the wind velocity, v w , set to 600 km/s 
and the mass-loading parameter (i.e., the ratio of the mass 
of gas given a velocity kick to that turned into newly formed 
star particles), rj, set to 2. This corresponds to using approxi- 
mately 40% of the total energy available from supernovae for 
a Chabrier (2003) IMF, which is assumed by default. This 
choice of parameters results in a good match to the peak 
of the star formation rate history of the universe (Schaye et 
al. 2009), but otherwise is somewhat arbitrary. However, the 
REF run is typical, in terms of the galaxy group properties 
it produces, of all of the OWLS runs that implement super- 
novae feedback only with efficiencies ^ 1 (i.e., varying rj and 
v w does not greatly affect the resulting group properties for 
constant ijv^,). This will be demonstrated in a forthcoming 
paper. 

The AGN run also includes kinetic supernovae feedback 
and uses the same choices for v w and r\. The prescription for 
the growth of BHs and feedback from AGN that we use, 
which is a substantially modified version of that introduced 
by SDH05, is described in detail in BS09. A basic description 
is as follows. Black holes grow by accretion of surrounding 
gas and by mergers with other black holes. In the simu- 
lations analysed in the present study, the Bondi radius is 
generally not resolved. This lack of resolution will result 
in an underestimate of the gas density near the BH, and 
therefore an underestimate of the BH accretion rate. The 
simulations particularly underestimate the Bondi rate if the 
gas is multiphase, since we impose an effective equation of 
state for gas with uh > n*. To compensate for this effect, 
the calculated Bondi-Hoyle-Lyttleton rate is multiplied by a 
coefficient a, as was done originally by SDH05. In the model 
of SDH05 (and the many studies that have used this model), 
a is fixed typically at a value of ~ 100 — 300, whereas BS09 
instead have a vary as a power law that depends on the 
local gas density. The power law index, f3, is fixed to 2 for 
the model we examine in the present study and the relation 
is normalised such that a — 1 when the local gas density 



is equal to the density threshold for star formation^. This 
approach has the same number of free parameters (one) as 
the implementation of SDH05, but has the advantage that 
the accretion rate does not exceed the Bondi rate when the 
Bondi radius is resolved in the simulation. In any circum- 
stance, the gas accretion rate may not exceed the Eddington 
rate, and most of the BH mass is assembled with accretion 
rates at or close to the Eddington rate. 

A certain fraction of the rest mass energy of the accreted 
gas, e, is used to heat a number of randomly selected neigh- 
bouring gas particles, 7iheat> by raising their temperatures 
by an amount ATh oa t. (Note that the BHs store feedback 
energy until it is sufficient to heat ^hcat by AT heat .) This 
efficiency, which is actually the product of two separate effi- 
ciencies, namely the radiative efficiency of the BH accretion 
disk and the fraction of the emitted energy that is assumed 
to couple to the gas, is a free parameter of the model. BS09 
find that a value of e = 0.015 yields a good match to the 
z = relations between BH mass and stellar mass and veloc- 
ity dispersion and the z = cosmic BH density. We therefore 
fix e at this value. 

Aside from the efficiency, the AGN run is characterised 
by the following 3 main parameters. The index in the power 
law scaling of the Bondi accretion coefficient, a, is set to 
/? = 2, rihcat = 1, and, ATh cat = 10 s K. In a forthcoming 
paper, we show the properties of galaxy groups simulated 
with AGN feedback are relatively insensitive to changes in 
/3 and rthcat , but that the results are sensitive to the temper- 
ature increase, particularly if it is comparable to, or smaller 
than, the virial temperature of the galaxy group, which for 
our sample is ~ 10 7 K. The reason for this is obvious, if the 
gas is not heated to a higher temperature than the surround- 
ing ambient medium, then its entropy and cooling time will 
remain low. This will rapidly lead to a pile-up of low en- 
tropy material in the central regions of groups and clusters, 
violating observational constraints. We have found that a 
temperature increase of ATh cat > 10 s K resolves this prob- 
lem on the group scale. 

It should be noted that, as individual gas particles are 
selected randomly for heating, the mode of feedback from 
BHs developed by BS09 will not necessarily produce large 
sets of cavities (or 'bubbles') of the kind frequently ob- 
served in massive nearby cool core clusters. (Thus, the model 
of BS09 does not distinguish between 'quasar' and 'radio' 
modes, unlike the model developed by Sijacki et al. 2007). 
We are currently developing models in which bi-polar bub- 
bles are inflated, similar to that of Sijacki et al. (2007), to see 
if this results in significant changes in the ability to match 
observations on the group scale. Qualitatively, though, we 
would expect the results to be similar, as to be successful 
in staving off catastrophic cooling in specific directions any 
such model will likely have to heat the gas over the full 4-zr 
steradians in a time-averaged sense (which happens by con- 
struction in the BS09 model of BH feedback). 



2 As discussed in BS09, any simulation that resolves the Jeans 
scales automatically resolved the Bondi radius of any BH whose 
mass exceeds the particle mass. Thus, normalising the power law 
between a and p in this way ensures that a goes to unity when 
the Bondi radius is resolved. 
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2.1.1 Analysis of galaxy groups 

For each z — galaxy group with M200 ^10 Mq , of which 
there are approximately 200 in both runs, we produce sim- 
ple mock observations. For the hot gas we generate synthetic 
X-ray observations, whereas for the stars we produce surface 
brightness maps in a variety of optical/near-infrared bands 
(e.g. Ks, B, V). These maps are used to i) identify the peak 
of the X-ray surface brightness distribution, which we desig- 
nate as the centre of the system when constructing 2D or 3D 
radial profiles (such as density, temperature, element abun- 
dances, etc.), and ii) for the stars, to identify the central 
brightest galaxy (the CBG) and its centre. 

In the case of hot gas properties, such as mean tem- 
perature or gas density profiles, we emission-weight the gas 
propertied so that an approximate like-with-like compari- 
son can be made directly with the observations. The X-ray 
luminosity of the j th gas particle is computed as: 



Lx,j = n e>j n H ,jAjVj 
X e (Z 3 ) 



(1) 



Pi 



v, 



X e (Zj 



Pj 



[X e (Zj) + XiiZ^f n(Z-)m H \i{Z-)m s 



A, 



where p is the particle gas density, m ga s is the particle mass, 
the volume is V — m gas /p, n e , «h and ni are the num- 
ber densities of electrons, hydrogen, and ions, respectively, 
X e = n e /nn, Xi = rn/nn, Z is the metallicity, p is the mean 
molecular weight, run is the mass of a hydrogen atom, and 
A is the cooling function in units of ergs cm 3 s _1 (integrated 
over some appropriate passband, such as 0.5-2.0 keV). We 
compute A by interpolating a pre-computed table generated 
using the Astrophysical Plasma Emission Codqj (APEC 
vl.3.1; see http://cxc.harvard.edu/atomdb/). APEC cool- 
ing rates are computed on an element-by-element basis and 
summed to yield the total cooling rate of each particle; i.e., 



(2) 



3 Physical quantities derived from real X-ray observations are 
necessarily emission- weighted (i.e., in an annulus the most lu- 
minous gas will contribute the most to the signal). This is why 
we emission-weight all of the quantities from the simulations. In 
practice, though, noticeable differences in the mass- and emission- 
weighted profiles only occur in the outskirts of the groups, which 
are typically not observable anyway. Note also that the emission- 
weighted temperature that we calculate is defined in precisely the 
same way as by Sun et al. (2009), to which many of our compar- 
isons are made. The emission-weighted temperature differs from 
the Mazzotta et al. (2004) definition of "spectroscopic" tempera- 
ture typically by only 10%. 

4 To maintain strict consistency with the implementation of ra- 
diative cooling in the simulations, it would be more appropriate 
to use cooling rates predicted by the CLOUDY software package 
(Ferland et al. 1998). The APEC, however, is more widely used 
in the analysis of X-ray data, which is why we have adopted it 
here. We have compared the rates predicted by the APEC against 
those predicted by the CLOUDY and MEKAL (another widely 
used plasma model) codes and, reassuringly, we find negligible 
(less than a few percent) differences between the predictions in 
the regime in which we are interested (i.e., hot, dense gas). 



where T is the gas temperature and \j,k{Tj) is the cooling- 
function for element species k for the j th particle. The sum- 
mation is performed over the 11 most important elements 
for cooling (H, He, C, Ca, N, O, Ne, Mg, S, Si, Fe), which are 
individually and self-consistently tracked during the simula- 
tion. 

Optical luminosities are computed by treating each star 
particle as a simple stellar population (SSP), which is rea- 
sonable as the particle mass is comparable to that of a star 
cluster. The simulations adopt a Chabrier IMF and store 
the age and metallicity of the particles. We use this informa- 
tion to compute a spectral energy distribution for each star 
particle using the GALAXEV model of Bruzual & Char- 
lot (2003). The luminosity is obtained by integrating the 
product of the SED with an appropriate transmission filter 
function (e.g., B, V, or K s filters). 



3 RESULTS 

Below we make detailed comparisons of the two OWLS runs 
with a wide variety of observational data. In Section 3.1 we 
examine the entropy and temperature distributions of the 
hot gas. The gas mass fractions of groups are presented in 
Section 3.2. In Section 3.3 we analyse the 'classical' X-ray 
mass-temperature and luminosity-temperature scaling rela- 
tions. In Section 3.4 we examine the overall star formation 
efficiency of galaxy groups, the mean ages of the CBGs, and 
the fraction of CBGs that are currently forming stars. Fi- 
nally, in Section 3.5 we analyse the gas- and stellar-phase 
metallicities of groups. 



3.1 Entropy and temperature distributions 

3.1.1 Entropy 

We begin by examining the entropy distribution of the gas 
in the REF and ^GA^OWLS runs and how it compares with 
that of observed galaxy groups. As is common in X-ray as- 
tronomy, we define the 'entropy', 5, as ksT/n^ 3 and use 
units of keV cm 2 . The observational entropy S is related to 
the thermodynamic specific entropy s via s cx In S and, like 
the thermodynamic entropy, S will be conserved in any adia- 
batic process. Therefore, unlike the temperature or density, 
which may be raised/lowered by any processes that com- 
press/expand the gas, the entropy will maintain a record of 
the thermodynamic history of the gas (see Voit et al. 2003; 
Voit 2005). 

In Fig. 1, we show the median 3D emission- weighted 
entropy profiles for the REF and AGN runs. This is the me- 
dian profile for all galaxy groups with masses in the range 
13.25 ^ log lo (M5oo/M0) ^ 14.25 [approximately 70 groups 
with a median mass of log 10 (M5oo/Mo) » 13.5], where M500 
is the total mass within rsoo, which is the radius that en- 
closes a mean density of 500 times the present-day criti- 
cal density of the universe. Shown for comparison are the 
hatched regions, which represent the 20th and 80th per- 
centiles of the XMM-Newton sample of Johnson, Ponman, & 
Finoguenov (2009) [approximately 20 groups with a median 
mass of log lo (M 5 oo/M0) « 13.5] and the Chandra sample 
of Sun et al. (2009) [approximately 40 groups with a me- 
dian mass of log 10 (M5oo/MQ) ~ 13.8]. Also shown (dotted 
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Figure 1. Radial entropy distribution of the gas in the REF 
and AGN OWLS runs. Shown are the median emission- weighted 
3D entropy profiles for the REF (solid blue curve) and AGN 
(dashed red curve) OWLS runs for all groups with 13.25 ^ 
log lo (M5oo/M0) ^ 14.25. The hatched black and green re- 
gions represent the XMM-Newton data of Johnson, Ponman, &: 
Finogucnov (2009) and the Chandra data of Sun et al. (2009), 
respectively. The hatched regions enclose the 20th and 80th per- 
centiles of the data. The simulated entropy profiles are remarkably 
similar in the central regions, implying that BH feedback, while 
preventing excessive star formation in the CBG (see Section 3.4), 
does not significantly disturb the hot gas. 



line) is the power law fit of Voit, Kay, & Bryan (2005) to 
the entropy profiles of a sample of groups and clusters sim- 
ulated with non-radiative physics (i.e., the self-similar an- 
swer). The radial coordinate has been normalised by rsoo. 
Typically, rsoo ~ 0.65r2oo- The entropy has been normalised 
by S500, the 'virial entropy', which is defined as: 



5*500 (z = 0) = 



"■e,500 



12/3 



(3) 



GM 5 oofirn H 



2r 500 [500f bPcr u{z = 0)/(fi e m H )]^ 3 



where ft = f2f,/f2 m = 0.1756, and fi e is the mean molecular 
weight per free electron. Note that, S500 is not the entropy 
at rsoo, even in the case of the self-similar model. This is be- 
cause S500 is defined in terms of the mean electron density 
within rsoo (as opposed to the density at rsoo) for a system 
with the universe baryon fraction, and the virial tempera- 
ture, fcsTsoo = G M500 fimH /2r 500, which is not equivalent 
to the gas temperature at rsoo for non-isothermal mass dis- 
tributions (for an NFW distribution, e.g., T(rsoo) ~ O.6T500; 
e.g., Loken et al. 2002). S500 therefore depends only on the 
halo mass, which is dominated by dark matter, and therefore 
is insensitive to the thermodynamic state of the gas. 

Both the observed and simulated groups exhibit a sig- 
nificant degree of 'excess entropy' with respect to the self- 
similar scaling, generally confirming the results of previous 
studies (e.g., Finoguenov et al. 2002). This is a tell-tale 



sign of the impact of non-gravitational physics (i.e., cooling 
and feedback) on the properties of the intragroup medium. 
Within rsoo (at least), all of the gas has been significantly af- 
fected by non-gravitational physics. This is contrary to what 
is seen in the most massive clusters (Tx > 4 keV), where at 
intermediate radii the gas reverts to the self-similar profile 
(McCarthy et al. 2008; Pratt et al. 2009). 

At small radii (r < O.lrsoo), the two OWLS runs yield 
profiles that are very similar and are consistent with the 
data, although the simulated gradient may be slightly too 
steep with respect to the data. The similarity of the profiles 
from the two runs at small radii is somewhat surprising, as 
in the AGN run the time since the last feedback outburst 
from the central BH is typically short (~ 10 7-8 yrs). This 
implies that, although the heating from the AGN is suffi- 
cient to prevent excessive star formation in the CBG (see 
Section 3.4), it does not significantly modify the emission- 
weighted properties of the hot intragroup gas at radii out 
to ~ O.lrsoo. This may be consistent with the observational 
findings of Jetha et al. (2007), who found very little differ- 
ence in the entropy profiles between groups with and without 
bright central radio sources. The similarity of the simulated 
group profiles may stem from the fact that the AGN heats 
only small packets of gas at a time, which are able to float 
buoyantly out to larger radii without doing much work on 
the surrounding gas (see discussion in McCarthy et al. 2008). 

At larger radii, however, the two models show an inter- 
esting departure from one another (with the AGN run hav- 
ing more excess entropy) that effectively brackets the obser- 
vational data. Note, however, that since the median masses 
and mass distributions of the simulated and observed group 
samples are not identical it is difficult to assess based purely 
on Fig. 1 which model gives the best match to the data. We 
address this mass dependence directly below. 

While enlightening, median profiles remove informa- 
tion about trends with system mass and scatter at fixed 
mass. One way to examine the trends with mass and scat- 
ter is to plot the entropy at several fixed characteristic 
radii versus system mass, or some proxy for system mass 
such as the mean X-ray temperature (see, e.g., Ponman 
et al. 2003). In Fig. 2, we show the entropy measured at 
0.15r5oo, r25oo, and rsoo vs. mean X-ray temperature (note 
that r2soo ~ 0.45rsoo) for all « 200 simulated groups with 
M200 > 10 13 Mq . We make a like- with- like comparison with 
the data of Sun et al. (2009) by emission-weighting the en- 
tropy of the simulated groups and computing 'cool core cor- 
rected' emission-weighted mean temperatures (see Section 
3.3 for a discussion of this 'cool core correction'), by excising 
the gas within projected radii of 0.15rsoo when computing 
the mean temperature. 

Consistent with the median profile results, we find that 
the two models are very similar at small radii (left panel) 
and consistent with the data. Interestingly, there is a much 
larger degree of scatter in the entropy at small radii in the 
AGN run, which will manifest itself later as increased scat- 
ter in the luminosity-temperature relation (see Section 3.3) 
and which perhaps also gives rise to scatter in the present- 
day star formation rate of the CBG (see Section 3.4). At 
larger radii, the two models begin to depart, with the AGN 
model yielding a marginally better match to the Sun et 
al. data, particularly at intermediate radii. Note that this 
slightly better match to the data is not as apparent in the 
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Figure 2. Entropy at several different characteristic radii as a function of cool core-corrected emission-weighted system temperature. 
Left: Entropy at 0.15rsoo- Middle: Entropy at T2500- Right: Entropy at rsoo- The error bars represent the measurements of Sun et al. 
(2009). The differences between the two OWLS runs are largest for systems with Tx j$ 1 keV, for which there are currently few entropy 
measurements. 



median profiles plotted in Fig. 1, which is a result of small 
differences in the median masses and mass distributions of 
the simulated and observed group samples that were stacked 
to create the median profiles. 

The difference between the AGN and REF runs at in- 
termediate/large radii is greatest for small system masses 
where presently there is not much observational data avail- 
able. Extending the observed samples to lower masses would 
therefore be very useful for helping to discriminate between 
the models, although the intrinsically low surface bright- 
nesses of such systems will not make such measurements 
easy by any means. 

Based on the above, we find that the effect of AGN feed- 
back on the entropy distribution of the intragroup medium 
is rather modest (even though, as we will see, AGN have 
a large effect on the ability of the gas to form stars), with 
the largest effects actually occurring at large radii. However, 
such similar behaviour in the entropy is actually what is re- 
quired to match the observations, as models with no AGN 
feedback (such as the REF model) match the observed en- 
tropy distribution relatively well. On the basis of entropy 
alone, therefore, there is not a compelling case for the need 
for AGN feedback. 

Our findings are broadly consistent with the results of 
several previous studies (e.g., Muanwong et al. 2002), who 
found that either excessive cooling or strong feedback can 
raise the entropy of the ICM to a level comparable to what 
is observed (the former mechanism does so by selectively re- 
moving the lowest-entropy gas via star formation; see Bryan 
2000; Voit & Bryan 2001). However, it is still somewhat of a 
surprise that the detailed radial entropy profiles of the REF 
and AGN runs are so similar. From these results alone it 
is not possible to distinguish what physical mechanism has 
caused the loss of the low entropy material. 

3.1.2 Temperature 

In Fig. 3, we show the median 2D (projected) emission- 
weighted temperature profiles for the REF and AGN runs. 
These are the median profiles for all galaxy groups with 
masses in the range 13.25 ^ log 10 (M5oo/Mo) ^ 14.25. 
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Figure 3. Median emission-weighted 2D (projected) tempera- 
ture profiles. The shaded black and green regions represent Chan- 
dra data of Rasmussen & Ponman (2009) and Sun et al. (2009), 
respectively, enclosing the 20 th and 80 th percentiles of their data. 
The simulated profiles depart within r < O.lrsoo due to the ex- 
cessive accumulation of cold baryons in the centres of groups in 
the run with supernova feedback only. 



Shown for comparison are the hatched regions, which rep- 
resent the 20th and 80th percentiles of the Chandra sample 
of Rasmussen & Ponman (2009) [15 groups with a median 
mass of log lo (M5oo/M0) ~ 13.5] and the Chandra sample 
of Sun et al. (2009) (see Section 3.1.1 for description of this 
sample). 

Through hydrostatic equilibrium (HSE), the gas tem- 
perature distribution can be thought of as a manifestation 
of the underlying entropy profile and the depth of the grav- 
itational potential well (see discussion in Voit 2005). If dark 
matter dominates the gravitational potential well and if its 
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Figure 4. The gas mass fraction of galaxy groups as a function of cool core-corrected emission-weighted system temperature. Left: The 
gas mass fraction within T2500 [i- e -> M gas (r25oo)/M250o] norma li se| i by the universal value, /[, = f2(,/f2 m = 0.18 vs. cool core-corrected 
emission- weighted temperature. Right: The analogous plot for r^oo - The hatched region represents the spread in gas mass fractions within 
r 500 f° r massive clusters (see McCarthy et al. 2007). The error bars in both panels represent the CTi-ancira-derived measurements of Sun 
et al. (2009). Energy input from supermassive BHs significantly reduces the gas mass fractions of groups with M500 j5 10 14 Mq, bringing 
them into agreement with the observations. 



distribution is relatively insensitive to what the baryons are 
doing, one expects that similar entropy distributions will 
give rise to similar temperature distributions. Indeed, we 
see that at large radii (r > O.lrsoo) the temperature distri- 
butions from the two runs are similar, and characterised by 
a mild negative gradient that is also seen in the observed 
groups. The ^GiVrun yields a slightly better fit to the data 
at large radii. 

Interestingly, the runs produce rather different temper- 
ature profiles within the central regions (r < O.lrsoo)- Since 
the entropy distributions at small radii are very similar (Fig. 
1) and since the gas is very close to HSE in both runs at 
small radii (we have verified this), this can only signal that 
the temperature difference is due to a difference in the depth 
of the potential wells in the very central regions. In partic- 
ular, increasing the amount of mass in the central regions 
will give rise to additional compressional adiabatic 'heating' 
(a rise in temperature) as gas flows into the centre. The 
higher central temperatures in the REF run therefore im- 
plies that there has been a much larger degree of accumula- 
tion of (cold) baryons (likely accompanied by a contraction 
of the dark matter halo) at the centres of groups in this run, 
than in the AGN run. 

The observed groups typically show evidence for a pos- 
itive temperature gradient in the very central regions. If one 
assumes that the dark matter distribution of the observed 
systems is similar to the simulated groups (e.g., Gastaldello 
et al. 2007), and given that the entropy distribution is sim- 
ilar to the simulated groups, the implication is that fewer 
cold baryons have accumulated in the central regions of real 
groups than predicted by the simulation that includes su- 
pernova feedback only. This is a prediction we explicitly 



test in Section 3.4. The run that includes feedback from 
BHs, however, yields central temperatures that are similar 
to observed, although the positive gradient appears to be 
somewhat too shallow. Based on this better agreement, we 
therefore predict that this run should also yield much more 
sensible stellar masses within the central regions of groups. 

Finally, based on the above we would argue that the 
presence/absence of a central positive gradient in the tem- 
perature profile is a poor proxy for the importance of ra- 
diative cooling, at least for simulations. The entropy would 
be a better quantity for assessing the importance of cooling, 
as it is unaffected by adiabatic compression. Which entropy 
threshold one adopts for characterising 'cool core' and 'non- 
cool core' systems is somewhat arbitrary, but this is unim- 
portant as long as the same threshold is adopted for both 
the simulations and observations. 



3.2 Gas mass fractions 

Another important property that is sensitive to the efficiency 
of radiative cooling and feedback processes is the gas mass 
fraction, / gas (i.e., the ratio of the mass of hot gas to to- 
tal mass within some characteristic radius). Until recently, 
it was only possible to directly measure the gas mass frac- 
tions in groups out to relatively small radii. However, deeper 
observations, particularly with Chandra, are now allowing 
measurements of the gas mass fractions of groups as far out 
as rsoo (e.g., Gastaldello et al. 2007; Rasmussen & Ponman 
2007; Sun et al. 2009). This gives a much larger lever arm 
for discriminating between thermal histories of the intra- 
group medium and it also allows one to compare groups to 
clusters at the characteristic radius where clusters show evi- 
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dence for convergence in / gas (see, e.g., Vikhlinin et al. 2006; 
McCarthy et al. 2007). 

In Fig. 4, we plot the gas mass fractions of groups in the 
REF and AGN OWLS runs within two characteristic radii, 
?"2500 and rsoo, vs. cool core-corrected emission-weighted 
temperature. For comparison we also show Chandra- derived 
gas mass fractions from Sun et al. (2009). The observed 
groups show very strong trends in / gas with system tem- 
perature (mass) at both r25oo and rsoo, with the lowest 
mass systems having the lowest gas mass fractions. The pre- 
dicted / gas — Tx relation from the A GN run is in excellent 
agreement with the data at both r2soo and rsoo- The REF 
(no AGN feedback) run, by contrast, shows a much milder 
trend in / gas with system temperature and it is only for the 
most massive systems that simulated groups resemble the 
observed ones. 

These results appear to agree very well with those de- 
rived recently by Puchwein et al. (2008) and Fabjan et al. 
(2009), who also implemented BH growth and AGN feed- 
back in self-consistent 'zoomed' cosmological simulations of 
a number of groups and massive clusters. This agreement 
comes in spite of the many differences in detailed imple- 
mentations of cooling, star formation, chemo dynamics, SNe 
winds, and AGN feedback in the OWLS simulations and 
those of Puchwein et al. (2008) and Fabjan et al. (2009). 

A simple qualitative explanation for these results is that 
additional energy input from supermassive BHs is sufficient 
to blow baryons out of low mass haloes or, alternatively, to 
raise the entropy of the gas to a high enough level to pre- 
vent it from ever being fully accreted by the system. The 
effect of feedback is largest, in a relative sense, for low mass 
haloes, which have the smallest gravitational binding ener- 
gies. The gas mass fractions of the two OWLS runs start to 
overlap at the highest system masses, for M500 > 10 14 M©. 
This therefore likely signals the mass scale where even su- 
permassive BHs do not inject sufficient energy into the gas 
to prevent systems from accreting nearljjf] their full com- 
pliment of baryons. This agrees remarkably well with the 
semi-analytic AGN model results of Bower et al. (2008) . 

3.3 Classical X-ray scaling relations 

3.3.1 Mass-temperature relation 

Fig. 5 shows the predicted relations between the total mass 
within rsoo and the cool core-corrected emission-weighted 
mean temperature for the REF and AGN inns. Also shown 
is the observed relationship of Sun et al. (2009) , derived by 
applying a hydrostatic analysis to Chandra data of a sample 
of approximately 40 groups. Overall, both models reproduce 
the observed relationship relatively well, especially the slope 
and scatter. These results are consistent with the findings of 
several previous studies based on cosmological simulations, 
including Borgani et al. (2004) and Nagai et al. (2007). 

There is perhaps a hint of a slight ~ 10% deviation in 
the normalisation of simulated and observed relations. One 
possible explanation for this deviation is that the observed 

5 To be exact, one must also account for the stellar mass fraction 
to find the halo mass at which systems become 'baryonically- 
closed'. This occurs at approximately M500 ~ 3 X 10 14 Mq in the 
AGN run (and for M500 > 10 13 M Q in the REF run). 
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Figure 5. The total mass within rsoo v s. cool core-corrected 
emission-weighted temperature. Error bars represent the Chandra 
measurements of Sun et al. (2009). The two OWLS runs predict 
very similar relations. 

masses may be biased slightly low due to a break down in 
the assumption of strict HSE. Nagai et al. (2007) have shown 
that if M500 is measured by applying the assumption of HSE 
to cosmologically-simulated groups, they retrieve a M500 — 
Tx relation that matches the observed relation of massive 
clusters better than in the case where they use the true 
masses from the simulations. These authors have shown that 
non-thermal pressure support, in the form of turbulent gas 
motions, is present at the ~ 10 — 20% (of thermal pressure) 
level at radii of rsoo and beyond and this is what leads to 
the slight bias in the observed mass estimates. 

It may at first sight be somewhat surprising that the 
AGN and REF runs have such similar mass-temperature re- 
lations. This is primarily a consequence of the fact that the 
models produce fairly similar entropy distributions and, per- 
haps more importantly, that the X-ray temperatures have 
been 'cool core-corrected'. Cool core correction is a stan- 
dard procedure designed to reduce the scatter in scaling 
relations between some observable and system mass, usu- 
ally for the purposes of yielding more precise constraints 
on cosmological parameters derived from tests involving the 
use of cluster masses (e.g., the cluster mass function). Typ- 
ically, the correction involves excluding the central regions 
(the cool core, if one is present) when computing the mean 
system temperatur^f). As we demonstrated in Section 3.1.2, 



6 The term 'cool core-correction' has also been applied to the 
procedure of calculating the luminosity after excising the emis- 
sion from the central regions or from integrating a parametric 
model (e.g., a beta model) that has been fit to the X-ray surface 
brightness at large radii (typically beyond the cooling radius) and 
extrapolated inward. However, when we apply it, we mean only 
the correction of the mean temperature, by computing over the 
projected radial range of 0.15 — l.Orsoo- 
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Figure 6. X-ray luminosity-temperature relation. Left: X-ray luminosity (0.5-2.0 keV band) vs. uncorrected emission-weighted temper- 
ature. Right: X-ray luminosity (0.5-2.0 keV band) vs. cool core-corrected emission-weighted temperature. Points with error bars in the 
left hand panel represent the ROSAT measurements of Helsdon & Ponman (2000) and Mulchaey et al. (2003) (uncorrected tempera- 
tures), while points with error bars in the right hand panel represent the ROSAT measurements of Osmond & Ponman (2004) (corrected 
temperatures). The AGN run yields a better match to both the corrected and uncorrected relationships, owing primarily to the reduced 
central gas densities in that run. 



the two models produce fairly different temperature profiles 
within the very central regions (r < O.lrsoo) of the groups. 
Uncorrected scaling relations, therefore, are more useful for 
probing the astrophysics in the cores of groups and clusters. 
This is illustrated more clearly below. 



3.3.2 Luminosity-temperature relation 

In Fig. 6 we plot the 0.5-2.0 keV X-ray luminosity vs. both 
the uncorrected (left panel) and cool core-corrected (right 
panel) emission-weighted temperatures. For comparison, we 
show in the left panel the ROSAT measurements of Hels- 
don & Ponman (2000) and Mulchaey et al. (2003) (temper- 
atures uncorrected) and in the right panel the R OS AT mea- 
surements of Osmond & Ponman (2004) (cool core-corrected 
temperatures) . 

Focusing first on the left panel of Fig. 6, we see that the 
uncorrected Lx — Tx relation from the A GN run is slightly 
steeper than that from the REF run. This is a consequence 
of the reduced gas densities in the former. Neither run, how- 
ever, produces a relation that is as steep as the observed one, 
although the AGN run clearly performs better than the run 
with supernova feedback alone. 

Interestingly, the scatter in the A GN run is larger than 
that in the REF run, with the former having a small popu- 
lation of outliers (with Tx > 1 keV and Lx % 10 42 ergs/s) 
that lie well off the mean relation. What is the origin of the 
increased scatter and the outliers in the A GN run? 

In the right hand panel of Fig. 6, we show the luminosity 
- cool core-corrected temperature relation. We see here that 
the scatter in the AGN run has been reduced in general 
and also that there are no significant outliers off the mean 



relation. This immediately indicates that both the increased 
scatter of the main population and outliers in the left panel 
are due to variations in the temperature of the gas within 
0.15r5oo, since the luminosity for the simulated groups is the 
same in both panels. That feedback from AGN can increase 
the scatter in the luminosity-temperature relation through 
its effect on the mean temperature appears to be consistent 
with the observational findings of Croston et al. (2005). 

The outliers in the left panel of Fig. 6 must have very 
hot, yet dense (i.e., dense enough to affect the mean sys- 
tem temperature) gas in their central regions. It is likely 
that these systems have undergone very vigorous heating by 
their central supermassive BHs recently. The temperature 
increase of gas heated by the central AGN is set to 10 8 K in 
the current run, which is « 8.7 keV. If a sufficient number of 
particles have been heated, this would naturally explain the 
presence of the outliers. However, the gas must have been 
heated quite recently (within the past 10 8 yrs or so), so that 
the heated gas has not had sufficient time to rise buoyantly 
out of the central regions. 

Observationally, high temperature outliers appear to be 
rare (see data in left hand panel), which may signal a prob- 
lem with the current model of AGN feedback. It is possible 
that this problem, if real, may be simply resolved by in- 
creasing the temperature to which the gas is heated by the 
BH (to, say, 10 9 K), so that the heated gas is so hot that 
it does not shine in soft X-rays, and therefore will not con- 
tribute strongly to the emission-weighted temperature (see 
discussion in McCarthy et al. 2008). Increasing the temper- 
ature to which the gas is heated should also result in longer 
periods of time between outbursts, which should further de- 
crease the fraction of outliers at any given time. The longer 
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time between outbursts is required as more gas must be ac- 
creted by the BH in order for it to have sufficient energy 
to heat the gas to the higher temperature. Alternatively, a 
change in the detailed implementation of the AGN feedback 
may also resolve this issue. For example, Sijacki et al. (2008) 
implemented a mode of feedback which inflates bubbles con- 
taining a non-thermal relativistic population of cosmic rays. 
It would be interesting to see whether or not such a model 
avoids producing outliers like the ones we observe in the 
present OWLS A GN run. 

A comparison to the observational data in the right 
panel of Fig. 6 shows that the AGN run matches the lumi- 
nosity - cool core-corrected temperature relation relatively 
well, certainly much better than the run with supernova 
feedback only. Since the cool core-corrected temperatures 
are similar for the two models at fixed system mass (Fig. 
5), the reason the AGN run matches the data must be that 
the X-ray luminosity has been reduced by a decrease in gas 
density owing to feedback from BHs. 

Puchwein et al. (2008) and Fabjan et al. (2009) also 
compared their zoomed cosmological simulations to the ob- 
served luminosity-temperature relation of groups and clus- 
ters and likewise concluded that the addition of energy in- 
put from supermassive BHs results in a better match to the 
data, particularly for lower mass systems. However, Fabjan 
et al. found this came at the expense of producing too high 
of a central entropy on the group scale. We see no sign of 
this problem for the AGN run. Improved agreement with 
the observed luminosity-temperature relation has also been 
found when BH feedback is incorporated into semi-analytic 
models (e.g., Bower et al. 2008) and hybrid methods (Short 
& Thomas 2009). 



3.4 Formation of Stars 

3.4-1 Group luminosity and luminosity of the CBG 

Galaxy groups and clusters are special objects in the uni- 
verse, as they are the only objects for which it is presently 
possible to directly measure the total mass and distribu- 
tion of both stars and hot gas (which dominate the baryonic 
mass budget) out to a large fraction of the virial radius. Over 
the past decade or so, cosmological hydrodynamic simula- 
tions have had varying degrees of success when it comes 
to matching the observed properties of the hot gas. Con- 
sistently, however, it has proved very difficult to construct 
physical models that are able to match the stellar proper- 
ties of groups and clusters. In particular, observed groups 
and clusters show only a modest star formation efficiency 
of ~ 10% (i.e., only 10% of the baryonic mass is in the 
form of stars; see, e.g., Balogh et al. 2008), whereas it is 
not uncommon for cosmological hydrodynamic simulations 
to yield much largeiQ efficiencies of 30%-40%. This general 



' We note here that an important factor in comparing differ- 
ent cosmological simulations is whether or not the simulations 
in question included radiative cooling losses due to metal lines. 
As we will show in a forthcoming paper, simply changing from 
primordial cooling rates to including metal-dependent radiative 
cooling as well, can increase the stellar fraction by a factor of 2 
or more. 



over-efficiency problem has been termed the 'cooling crisis' 
of cosmological simulations (Balogh et al. 2001). 

It is important to recognise that this cooling crisis ap- 
peared in, and still only strictly applies to, cosmological sim- 
ulations that invoke feedback from supernovae winds (or no 
feedback at all) but not from supermassive BHs. The recog- 
nition in recent years that heat input from supermassive BHs 
can significantly affect the properties of the gas, and there- 
fore its ability to form stars, motivates us to re-examine the 
star formation efficiency of groups in cosmological simula- 
tions. 

In Figure 7 we plot the K-band stellar luminosity within 
f500 (left panel) and of the central brightest galaxy (CBG) 
(right panel) vs. the uncorrected emission-weighted mean 
X-ray temperature. Because it is in the near-infrared, the 
K-band luminosity is sensitive mainly to older stellar pop- 
ulations and therefore should be a much better tracer of 
stellar mass than, say, optical luminosities in the B or V 
bands. Furthermore, it should be relatively unaffected by 
dust extinction. Shown for comparison is the observational 
data of Lin & Mohr (2004; 2MASS K-band luminosities) and 
Horner (2001; ASCA uncorrected X-ray temperatures). The 
K-band luminosity of each star particle in the simulations is 
calculated using their formation times and metallicities and 
applying the simple stellar population models of Bruzual & 
Chariot (2003), with no dust correction. The CBG is identi- 
fied by centroiding on mock K-band surface brightness maps 
of the simulated groups. When quoting properties of the 
CBG below, we are referring to properties of the stars within 
30 hT x kpc of the CBG centre. For comparison, the typical 
effective radius of an observed CBG is 10-15 kpc (e.g., Shen 
et al. 2003). 

From Fig. 7 it is immediately apparent that the REF 
model produces far too many stars with respect to the ob- 
servations. There are too many stars within r-500, by approx- 
imately a factor of 4, and too many within the CBG, by a 
similar factor. This is the well known overcooling problem. 
Inclusion of feedback from supermassive BHs, however, has 
a dramatic effect on the ability of the gas to form stars. So 
much so that the star formation efficiency is lowered by al- 
most exactly the factor that is required to match the data. 
Note that the parameters of the AGN feedback were not 
tuned in any way to get this result (they were tuned to 
match the normalisation of the observed BH scaling rela- 
tions; BS09). Impressively, the slope of the observed relation 
is also reproduced relatively well by the simulations. Note 
that the outliers in the Lk—Tx relation of the AGN run are, 
as in Section 3.3.2, those systems which have anomalous cen- 
tral temperatures due to a very recent outburst - the scatter 
is much reduced when plotting vs. cool core-corrected tem- 
peratures, although no such temperatures are available for 
the Lin & Mohr K-band CBG sample. 

We point out that the fact that the run that incorpo- 
rates energy input from supermassive BHs results in CBGs 
with reasonable luminosities, while the run that does not 
over-predicts the stellar mass of the CBG, agrees well with 
our expectations based on the analysis of the temperature 
profiles in Section 3.1.2. In particular, the temperature pro- 
file rises all the way into the centers of groups in the run 
that neglects AGN feedback, indicating that the central po- 
tential wells are much deeper in that run. Therefore, X-ray 
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Figure 7. K-band luminosity within rsoo (Left) and of the CBG (Right) vs. uncorrected emission-weighted temperature. The points 
with error bars on both luminosity and temperature represent the data of Lin & Mohr (2004; K-band luminosity) and Horner (2001; 
X-ray temperature). The points with error bars only on the temperature represent the measurements of Rasmussen & Ponman (2009). 
Feedback from supermassive BHs reduces the efficiency of star formation by approximately a factor of 4, bringing the predicted Ljf — Tx 
relation into good agreement with the observations. 



observations offer a probe independent of optical methods 
of the efficiency of star formation in groups and clusters. 

As a word of caution, the overlap (in terms of system 
mass) between the observations and simulations needs to 
be increased in order to be more definitive about our con- 
clusions based on Fig. 7. Observationally, this is difficult 
as obtaining temperature measurements for the hot gas in 
very low mass groups will likely require very long integra- 
tion times and a very careful accounting of foreground and 
background X-ray sources. The accuracy in the assignment 
of group galaxy membership also becomes much more im- 
portant in the limit of low mass systems with few galax- 
ies. Alternatively, one can run simulations in larger periodic 
boxes to obtain more massive systems, for which observa- 
tions already exist. This has the disadvantage of focusing 
on relatively rare systems, so there is a worry about how 
general the conclusions obtained from comparing very mas- 
sive simulated and observed systems are. 

In addition, it should also be kept in mind that the lumi- 
nosities reported for the simulations represent the summa- 
tion of all star particles within rsoo , whereas the reported ob- 
served luminosities represent the summation of light above 
some surface brightness limit. It is therefore possible that 
the present AGN run may have even yielded a slightly too 
low star formation efficiency. There is currently a debate in 
the literature over what fraction of stellar mass/luminosity 
is in a diffuse, extended component (the intracluster light, or 
ICL for short), with estimates ranging from 10%-50% of the 
total light being in the ICL (Lin & Mohr 2004; Murante et al. 
2004; Zibetti 2005; Gonzalez et al. 2007; Balogh et al. 2008; 
McGee & Balogh 2009). The agreement with the central 
temperature profiles of observed groups (see Fig. 3) likely 
indicates that the fraction of stellar material in the central 



regions that is below the 2MASS surface brightness limit is 
small. In any case, it is very clear that there is insufficient 
mass/light in the ICL to reconcile supernovae feedback-only 
simulations with observations. More importantly, we have 
demonstrated that the inclusion of a physically-motivated 
model for feedback from supermassive BHs can significantly 
reduce the efficiency of cooling of gas destined to end up in 
groups (or group galaxies). This is consistent with findings 
of several previous studies, including Sijacki et al. (2007), 
Fabjan et al. (2009), and Puchwein et al. (2008, 2010). 



3.4.2 Mean age of the CBG 

The simulations should not only match the total present-day 
stellar mass, but they must also reproduce the star forma- 
tion histories of these systems, including their present-day 
star formation rates. Historically, supernovae feedback-only 
simulations have not only failed to yield reasonable stellar 
luminosity and mass functions, they have also failed to re- 
produce the colours of observed CBGs, with the simulated 
galaxies being too blue with respect to the observed sys- 
tems. This is an indication that recent star formation is 
too efficient in these simulations. Another way of expressing 
the same thing is that the mean stellar age of the CBG is 
too young compared to the observed ages. We now examine 
whether including feedback from supermassive BHs affects 
the mean stellar age of the CBG. 

In Fig. 8 we show histograms of the B-band emission- 
weighted (left panel) and mass-weighted (right panel) CBG 
age for the REF and AGN OWLS runs. Also shown, in the 
left panel, is the data of Loubser et al. (2009). These authors 
derived the stellar ages of approximately 50 CBGs in a sam- 
ple of nearby groups and clusters from long slit spectra using 
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Figure 8. B-band emission-weighted (Left) and mass-weighted (Right) stellar age of the CBG. Solid blue and dashed red histograms 
represent the REF and AGN OWLS runs, respectively. In the left panel, the dotted black histogram represents observational data of 
Loubser et al. (2009). Energy input from supermassive BHs significantly reduces the star formation rates of CBGs at late times, yielding 
much higher (older) mean CBG ages in good agreement with observations. 



the stellar population synthesis models of Maraston et al. 
(2003, 2004) with the Lick system of absorption line indices. 
(The advantage of using the absorption line indices is that 
they should not be affected by dust.) The derived quanti- 
ties are obviously emission-weighted and because they are in 
the optical a comparison to the B-band weighted simulation 
ages is appropriate (we have also compared to the V-band 
weighted simulation ages and find very similar results). 

Focusing first on a comparison between the models, we 
see that when the stellar ages are luminosity weighted there 
is a large difference in the mean age of the CBG between the 
two models, with the peak of the distribution being between 
^2-3 Gyr for the REF run and between tti 8 — 9 Gyr 
for the AGN run. The mass- weighted distributions are more 
similar and indicate a generally old population of stars in 
both runs. This is expected since most of the stellar mass 
forms at z ~ 1 — 3 (corresponding to the peak of the star 
formation rate density of the universe), even in simulations 
with inefficient feedback such as the REF run. These results 
indicate that the B-band luminosity weighted ages are very 
sensitive to recent episodes of star formation and that the 
star formation rate of CBGs in the REF run are presently (or 
have very recently been) much higher than those of CBGs 
in the simulation that includes AGN feedback. 

A comparison to the data of Loubser et al. (2009) con- 
firms our expectation that the REF model yields mean stel- 
lar ages that are far too young compared to real CBGs. 
While we should expect an improvement by including feed- 
back from supermassive BHs, agreement with the data is 
not necessarily automatic. It will obviously depend on when 
and where (that is, in what progenitors) the BHs heated the 
gas, and one can envisage not unreasonable scenarios where 
the integral of the star formation rate history is reproduced 
but not the detailed history itself. In spite of this, the AGN 



run reproduces the distribution of stellar ages of observed 
CBGs well, especially the peak of the distribution. The large 
scatter is also retrieved relatively well. 

On a qualitative level, our results confirm the findings 
of De Lucia et al. (2006) , who used a semi-analytic model of 
galaxy formation to show that the inclusion of a prescription 
of AGN feedback (see Croton et al. 2006) results in much 
older formation ages for massive elliptical galaxies. 

3.4.3 Present-day star formation rate of the CBG 

We now examine the present-day star forming properties of 
the CBGs in the two runs. It is important to note that it is 
entirely plausible that a model could match the K-band lu- 
minosity of the CBG without matching the present-day star 
formation rate. This is because most of the stellar mass was 
not formed recently out of a cooling flow, even in models 
where there is no feedback at all and cooling flows oper- 
ate uninhibited at the present day. Most of the stellar mass 
in groups and clusters was formed in galaxies at high red- 
shift (z ~ 1 — 3). Therefore shutting down a cooling flow 
today will not help with the overcooling problem of cosmo- 
logical simulations if excessive cooling is not also shut down 
in galaxies at high redshift. By the same token, getting the 
correct star formation rate of galaxies at high redshift does 
not guarantee prevention of cooling flows today. 

We show in Fig. 9 the present-day star formation rate 
of gas in the CBG vs. uncorrected X-ray temperature. A 
tight trend between these quantities exists for the REF run, 
with the SFR of the CBG increasing with system temper- 
ature (mass). At a temperature of 1 keV (M500 ~ 3 x 10 13 
Mq), CBGs in the REF groups are forming stars at a rate 
of ~ 50 — 100 Mq yr _1 . By contrast, when feedback from 
supermassive BHs is included, the star formation rates are 
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Figure 9. Present-day (z = 0) star formation rate of the cen- 
tral brightest galaxy (CBG) vs. uncorrected emission-weighted 
temperature. The red arrows denote upper limits on the star for- 
mation rate for CBGs in the AGN run with no star forming gas 
particles. The upper limit is the minimum resolvable SFR in the 
simulations, corresponding to the SFR of a single particle with 
a density equal to the threshold for star formation. AGN feed- 
back dramatically reduces the present-day star formation rates of 
CBGs, shutting off star formation altogether in most systems. 



drastically reduced and in most cases star formation is shut 
off altogether. 

What are the present-day star formation rates of real 
CBGs? This is actually not a straightforward question to an- 
swer. One does not directly measure the star formation rate, 
but instead uses observables (such as infrared, H-alpha, and 
radio luminosities) as proxies for the star formation rate. 
However, these observables do not correlate perfectly with 
the star formation rate and systematics such as dust or aper- 
ture corrections can be very important Also, one must take 
care to account for differences in the assumed stellar IMF. 
Thus, accurate derivation of absolute star formation rates 
for CBGs is non-trivial. One might expect, though, that the 
observations can be used to robustly determine the fraction 
of CBGs that are forming stars at a reasonable (detectable) 
rate. 

The fraction of CBGs forming stars has been shown 
to depend somewhat on sample selection criteria, as there 
is a known correlation between the cooling rate of the hot 
gas and whether or not there is (or has recently been) star 
formation in the central object (e.g., Crawford et al. 1999; 
Edwards et al. 2007; Bildfell et al. 2008; Rafferty et al. 2008; 
Wang et al. 2009). Thus, a sample containing a large frac- 
tion of cool core groups and clusters (as one might expect in 
the case of X-ray flux-limited samples) should be expected 



8 Hicks et al. (2009), for example, find very large discrepancies 
(up to an order of magnitude) between UV and IR-inferred SFRs 
for CBGs in a sample of 16 cool core clusters. 



to host a larger fraction of star forming CBGs than a sample 
that contains a small fraction of such systems. For example, 
Crawford et al. (1999) found that for a sample of 256 domi- 
nant galaxies in 215 nearby X-ray-selected groups and clus- 
ters, approximately 27 per cent have emission line spectra. 
More recently, Edwards et al. (2007) examined the optical 
spectra of CBGs in two separate samples; an X-ray-selected 
sample based on the NOAO Fundamental Plane Survey and 
an optically-selected sample based on the SDSS DR3 sur- 
vey. Edwards et al. found that only ~ 13% of the optically- 
selected CBGs contained line emission, whereas ~ 20% of 
the X-ray-selected CBGs contained detectable line emission. 
A similar trend, based on larger optical and X-ray samples, 
has been reported recently by Wang et al. (2009). While Ed- 
wards et al. and Wang et al. did not quantify what fraction of 
their optically-selected groups/clusters contained cool cores 
(as many of their optical systems have no existing pointed 
X-ray observations), there is no obvious reason why these 
systems should be biased for or biased against in an optical 
selection. 

The resulting fraction of CBGs that have robust line 
emission therefore likely lies between 10% — 30%. Assuming 
this line emission is due to ongoing star formation (as op- 
posed to, e.g., AGN), would imply that this fraction of CBGs 
are forming stars at a rate of at least a few solar masses 
per year, assuming a Chabrier IMF (M. Balogh, priv. com- 
munication). However, this likely represents an upper limit, 
as some of the optical line emission will not have originated 
from recent star formation. For example, O'Dea et al. (2008) 
examined a sample of 62 CBGs with optical line emission 
in the mid-infrared with the Spitzer Space Telescope and 
found that only half of the systems showed evidence for a 
detectable excess of mid-IR emission. Thus, the fraction of 
CBGs with ongoing star formation could actually be closer 
to only 5% - 15%. By contrast, 100% of the CBGs in the 
OWLS REF run are forming stars at a rate that would be 
detectable with present observations. When feedback from 
AGN is included, we find the fraction of CBGs that are form- 
ing stars at a rate of > 3 Mq yr _1 is reduced to only 14%, 
which is comparable to the observed fraction. 

Interestingly, for those systems that are forming stars in 
the AGN run, there is a large degree of scatter in the SFR. 
This may be linked with the large degree of scatter seen 
in the mean ages of the CBGs (Fig. 8) and in the central 
entropy of the hot gas (left panel of Fig. 2). In a future paper 
we will examine the physical origin of the scatter in the 
star forming properties of CBGs in simulations with AGN 
feedback. 



3.5 Metallicity 

The metallicity of gas plays a crucial factor in its ability 
to cool and form stars. The cooling rate obviously also has 
implications for the thermodynamic properties of the hot 
gas that has not yet turned into stars, particularly on the 
galaxy group scale where a large fraction of the X-ray emis- 
sion originates from line emission. In order to match the 
thermodynamic properties of the hot gas in groups in de- 
tail with a physical model, one should therefore also match 
the metallicity distribution of that gas. It is therefore im- 
portant that simulations allow for metal-dependent cooling 
losses and obviously for the physical mechanisms that are 
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Figure 10. Median emission-weighted 3D iron (Left) and silicon (Right) abundance profiles for the REF and AGN OWLS runs. Error 
bars represent the binned averages of the Chandra data of Rasmussen & Ponman (2009). Both the observational and theoretical radial 
profiles are normalised using the solar abundances of Asplund et al. (2005). Both models reproduce the iron data relatively well. The 
shape of the silicon profile agrees better with the run including AGN, but the predicted abundance may be slightly too high. 



responsible for producing and getting metals out into the 
gas. 

As described in Wiersma et al. (2009a), the OWLS simu- 
lations compute cooling rates on an element-by-element ba- 
sis for the 11 most important species for cooling. The timed 
release of these elements (into the gas surrounding star par- 
ticles) by Type II and Type la SNe, as well as from interme- 
diate mass stars ('AGB' stars), is tracked during the simula- 
tions (Wiersma et al. 2009b). Whether or not these metals 
will then get distributed through the intragroup medium 
will depend on whether the metal-enriched gas cools quickly 
and gets locked back up in stars or, if not, if feedback and/or 
dynamical processes associated with the group environment 
are sufficient to move metals from haloes of individual galax- 
ies into the intragroup medium. Thus, the gas- and stellar- 
phase metallicities of groups provide an important diagnos- 
tic of the efficiency of star formation, feedback, and dynam- 
ical processes (such as ram pressure and tidal stripping) in 
groups. We now examine how well the OWLS REF and AGN 
simulations are able to reproduce the observed metal content 
of groups. 

Before proceeding, it is important to note that nucle- 
osynthesis yields and SN la rates are both uncertain at 
the factor of ~ 2 level, even for a fixed stellar initial mass 
function (see Wiersma et al. 2009b). Therefore, at least at 
present, we argue that more weight should be placed on the 
successes/failures of simulations in reproducing the thermo- 
dynamic state of the gas and the overall baryonic mass frac- 
tions (discussed in Sections 3.1-3.4). Metal abundances can 
only be expected to match the observations to within a fac- 
tor of ~ 2. 



3.5.1 Gas-phase metallicity 

In Fig. 10 we plot the median emission-weighted iron and 
silicon profiles for the REF and AGN runs. For comparison, 
we plot the recent CTiandra-derived results of Rasmussen & 
Ponman (2009) for a sample of 15 bright galaxy groups. 

Focusing first on the profiles, we see that both OWLS 
runs produce similar large-scale metallicity gradients, with 
the iron (silicon) abundance reaching solar values (several 
times solar) at the centre and declining to a tenth (a few 
tenths) solar by rsoo- In fact, large-scale metallicity gradients 
such as these are a generic feature of all of the OWLS runs. 
In addition, we note that the profiles presented here appear 
to be very similar in shape to those found in the zoomed 
cosmological simulations of massive clusters in Tornatore et 
al. (2007) and Borgani et al. (2008). 

The shape and normalisation of the simulated iron pro- 
files are remarkably similar to the binned average profile 
found by Rasmussen & Ponman (2009). In terms of cooling 
of hot gas, iron is the most important species. It is there- 
fore fortunate that the simulations reproduce the gas-phase 
iron abundance, given the uncertainty in the nucleosynthesis 
yields and SN la rates. Both simulations appear to produce 
groups that are over-abundant in silicon. The shape of the 
observed silicon profile, at least in the central regions, is 
more similar to that of the AGN run, so an adjustment of 
the assumed Si yields would improve the fit to the data. One 
might naively expect that adjusting the yields could have 
important consequences for the thermodynamic state of the 
gas, and this would indeed be the case for iron. However, 
silicon is a far less efficient coolant than iron for the temper- 
ature range relevant for this study (see Fig. 6 of Wiersma et 
al. 2009a). Therefore, a change in the Si abundance would 
likely only have a very minor effect on the resulting gas 
properties. 
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Figure 11. Total mass of hot gas-phase iron as a function of 
cool core-corrected emission-weighted X-ray temperature. Error 
bars represent the measurements of Rasmussen & Ponman (2009). 
Both models agree with the data. 



The result that including feedback from supermassive 
BHs gives rise to ICM metallicity profiles whose shape is 
in better agreement with observed profiles, at least in the 
central regions, was also obtained recently by Fabjan et al. 
(2009). On the other hand, neither the REF nor the AGN 
models reproduce the shape of the observed silicon profile at 
large radii. Like Fabjan et al. we find that the Si/Fe ratio is 
nearly flat as a function of radius from roughly 0.2rsoo out to 
rsoo, whereas observationally Rasmussen & Ponman (2009) 
find a steep rise in Si/Fe beyond about 0.2rsoo, which is due 
to the flattening of the Si profile at large radii. However, 
observationally-inferred silicon abundances in the outer re- 
gions of groups are still controversial. As noted by Fabjan et 
al. (2009), several other recent observational studies based 
on Suzaku data (e.g., Sato et al. 2008, 2009) find no evidence 
for such a steep rise in Si/Fe at large radii (although note 
that typically these studies have been limited to « 0.3rsoo)- 
If the results of Rasmussen & Ponman (2009) are confirmed, 
this will signal a problem with the current models of AGN 
feedback that may not be rectified by a simple change of the 
adopted yields and/or SN la rates. 

Instead of stacking the groups to obtain a single median 
profile, we plot in Fig. 11 the total mass of iron within rsoo 
vs. system temperature (mass). When computing the total 
metal mass for the simulated groups, we sum only metals 
attached to gas particles that are hot, > 10 5 K, and there- 
fore X-ray-emitting. In accordance with the profiles results, 
the total mass of iron in the hot gas phase is in excellent 
agreement with the observations. We note that we have also 
computed the total mass of silicon within rsoo and compared 
it with the measurements of Rasmussen & Ponman (2009) 
and find similarly good agreement. This may seem surpris- 
ing, given the profiles plotted in Fig. 10, but note that most 
of the metal mass is located at large radii, beyond ~ 0.4rsoo, 
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Figure 12. B-band emission-weighted stellar-phase metallicity 
of the CBG, presented in logarithmic units and relative to the 
Sun. Solid blue and dashed red histograms represent the REF 
and AGN OWLS runs, respectively. The dotted black histogram 
represents observational data of Loubser et al. (2009). We adopt 
the solar abundances of Asplund et al. (2005) in normalising the 
simulations. 



where the AGN run reproduces the observed Si abundance 
relatively well. However, this agreement may be misleading, 
as the measured Si abundances are most uncertain at large 
radii and, as already noted, the absolute abundances in the 
simulations are uncertain at a factor of 2 level owing to un- 
certainties in the adopted nucleosynthesis yields and SNIa 
rates. 



3.5.2 Stellar metallicity 

It is a curious fact that the REF and AGN runs produce 
such similar gas phase metal masses (Figs. 10-11), given 
that they have very different total stellar masses (Fig. 7). 
This may seem to be a contradiction at first, but in fact the 
explanation for this behaviour is simple: the mass of metals 
locked up in stars is very different in the two runs. In Fig. 12 
we show the distribution of B-band emission-weighted CBG 
metallicities for the two OWLS runs. Also shown for com- 
parison are the observational determinations of Loubser et 
al. (2009) for a sample of approximately 50 nearby CBGs in 
groups and clusters. 

A comparison between the two simulations shows nearly 
an order of magnitude difference in the peaks of the two dis- 
tributions, with the CBGs in REF run being far more metal 
rich than those in the AGN run. The effect of feedback from 
supermassive BHs is therefore not only to quench the overall 
star formation rate, but also to prevent the metals that are 
produced from being locked back up in future generations 
of stars. The BHs heat local gas causing it to expand out of 
the dense, star forming-regions and thereby raise the frac- 
tion of the total metal mass that is in the gaseous phase. The 
efficient ejection of metals from star forming-regions is the 
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mechanism by which the AGN run yields similar gas- phase 
abundances to the REF run, in spite of it having a stellar 
mass fraction that is « 4 times smaller than the latter run. 
The BH heating also leads to increased mixing and a more 
uniform metal distribution. In this respect, our results agree, 
at least qualitatively, with those of Sijacki et al. (2007) and 
Fabjan et al. (2009). 

Loubser et al. (2009) have recently measured the 
emission-weighted stellar metallicity of a large sample of 
CBGs. It should be noted that the derived abundances are 
not fully self-consistent in terms of the assumed solar abun- 
dance, as the isochrones adopted by the Maraston et al. 
models (from Cassisi et al. 1998) assume Zq — 0.02, whereas 
the stellar library is based on the solar neighbourhood and 
the metallicities of individual stars are given relative to the 
Sun (P. Sanchez-Blazquez, priv. communication). For the 
simulated CBGs we assume the up-to-date solar abundances 
of Asplund et al. (2005), which gives Z Q = 0.012. 

A comparison of the simulated CBGs to the data of 
Loubser et al. (2009) reveals that neither run performs par- 
ticularly well compared to the data, with the REF run pro- 
ducing too metal rich and the AGN run producing too metal 
poor stars. If we were to adopt, instead of the Asplund et 
al. (2005) abundances, a solar metallicity of Zq — 0.02, this 
would shift the models to lower metallicities (in solar units), 
making the discrepancy larger for the ^GWrun and smaller 
for the REF. However, even with this shift the REF model 
does not provide a good match to the data. We have also 
compared the simulations to the [Fe/H] and [Mg/H] mea- 
surements of Carretero et al. (2007) for a sample of 27 el- 
liptical galaxies in 4 nearby massive clusters and found very 
similar results. 

The apparent discrepancy may indicate that the 
adopted yields are not correct in detail and/or that the 
metal ejection resulting from supermassive BH feedback in 
the current simulation may be too efficient. It is presently 
unclear what the specific cause of the discrepancy is. If the 
problem is related to over-efficient metal ejection, we expect 
that rectifying this problem would not significantly affect 
the excellent match to the gas phase metallicity, as long as 
the match to the observed stellar content is maintained. The 
reason for this is that only a small fraction of the baryons 
in observed groups are in the form of stars. Thus, only a 
relatively small amount of metals would need to be removed 
from the intragroup medium to raise the stellar metallicity 
to the observed level. 



4 DISCUSSION & CONCLUSIONS 

We have analysed two high resolution cosmological hydro- 
dynamic simulations (the REF and AGN runs) from the 
OWLS project. While both include galactic winds driven by 
supernovae, only one of the models includes feedback from 
accreting BHs (the AGN run). These runs are representa- 
tive of the OWLS runs with supernovae feedback only and 
runs with AGN feedback, respectively, in terms of the galaxy 
groups properties they produce (Paper III). We compared 
the runs to a wide range of observations in order to ascertain 
whether or not including BH feedback in cosmological sim- 
ulations yields a more realistic population of galaxy groups. 
This comparison has yielded the following findings: 



• The two models yield similar median entropy distribu- 
tions for the gas and reproduce the 'excess entropy' of ob- 
served groups with respect to the self-similar scaling. How- 
ever, there is a much larger degree of scatter in the central 
entropies of groups in the run that includes BH feedback. 
This, in turn, gives rise to a larger scatter in the luminosity- 
temperature relation of this model. 

• The two models produce similar temperature profiles 
beyond O.lrsoo, but within this radius the temperature of 
groups in the REF model continue to rise, in discord with the 
observations. This is due to the excessive accumulation of 
cold baryons in the centres of groups in the absence of AGN 
feedback, which significantly deepens the central potential 
well leading to compressional heating of the gas. 

• Supernovae-driven winds alone do not inject sufficient 
energy to reproduce the observed low gas mass fractions 
of galaxy groups. The additional energy input from BHs is, 
however, sufficient to reduce the gas mass fractions of galaxy 
groups to the observed level and also reproduces the steep 
increase in the gas fraction with system temperature. 

• Both models match the observed scaling between mass 
and cool-corrected temperature relatively well. This is be- 
cause beyond the central regions the gravitational potential 
wells are similar in both models, owing to the dominance of 
dark matter at large radii. 

• The reduced gas density in the groups that have been 
subjected to BH feedback results in lower X-ray luminosities 
at fixed halo mass/mean temperature. The resulting lumi- 
nosity - cool core-corrected temperature relation is in good 
agreement with observations. When the temperatures are 
not cool core- corrected, the scatter in the L-T relation is in- 
creased significantly for the AGN run, with several outliers 
lying well off the mean relation. This is due to recent large 
episodes of BH heating. Such outliers appear to be very rare 
in nature. We speculate that if the heating temperature were 
increased to a higher level than 10 s K the scatter would be 
reduced, as the heated particles would cease to emit most 
of their radiation in soft X-rays. Also, the length of time 
between outbursts would likely increase, as more mass will 
have to be accreted by the black hole in order to have suffi- 
cient energy to heat the gas to a higher temperature. 

• The AGN run yields stellar K-band luminosities of both 
overall group and the CBG that are in good agreement with 
observations. The inclusion of feedback from BHs also yields 
the correct emission-weighted mean age of CBGs (as well as 
the scatter in age) and the fraction of CBGs that are forming 
stars at the present day (~ 15%). In contrast, the simulation 
with supernova feedback alone produces CBGs that are far 
too massive and young (with very little scatter in mean age) 
and all of the CBGs in this run have significant present-day 
star formation rates. 

• Both models yield total hot gas-phase iron masses that 
agree with observed masses and both yield large-scale gra- 
dients that are similar to the observed profiles. However, 
while the simulations reproduce the observed slope and nor- 
malisation of the iron abundance profile, both appear to be 
over-abundant in silicon at small radii. The discrepancy is, 
however, not greater than a factor of two and therefore lies 
within the current uncertainty in nucleosynthesis yields, so 
at present this is not a significant problem for the simula- 
tions. 
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• Neither simulation yields the correct metallicity of 
stars, with the simulation with BH feedback producing too 
metal poor stars while the simulation without BH feedback 
produces too metal rich stars. This may signal that the 
present implementation of AGN feedback is too efficient in 
expelling metals from star forming regions and/or that the 
adopted yields are incorrect in detail. If it is the former, we 
note that since the mass in metals in the hot gas phase is 
much larger than that locked up in stars, only a relatively 
small fraction of metals would be needed to be removed from 
the intragroup medium to resolve this problem. 

Based on the above, the inclusion of feedback from su- 
permassive black holes significantly improves the ability of 
our cosmological hydrodynamical simulations to yield a real- 
istic population of galaxy groups. This agrees qualitatively 
with the findings of Sijacki et al. (2007), Puchwein et al. 
(2008, 2010), and Fabjan et al. (2009), who implemented 
BH growth and feedback in cosmological zoomed simula- 
tions. We note, however, that our conclusions are based on 
a comparison of a much larger sample of simulated groups 
to a wider range of observations. In addition, the OWLS im- 
plementation of BH feedback appears to have yielded bet- 
ter matches to several different key observables than those 
found in some of these previous studies. For example, Fabjan 
et al. (2009) and Puchwein et al. (2008, 2010) report stellar 
mass fractions within rsoo that are a factor of ~ 2 higher 
than observed and what we obtain from the AGN model. 
In addition, Fabjan et al. (2009) find their galaxy groups 
have too high entropies at r25oo- Given the many differences 
in the detailed implementations of cooling, star formation, 
chemodynamics, and feedback in the OWLS simulations and 
those which are based upon the model of SDH05, it is not 
trivial to determine what physics is most responsible for the 
better match to the data. One possibility is the simulations 
of Sijacki et al. (2007) and Puchwein et al. (2008) did not in- 
clude metal-line cooling which is expected to be important. 

We have highlighted a number of areas where the cur- 
rent model of AGN feedback could be improved, includ- 
ing reproducing the scatter in the (uncorrected) luminosity- 
temperature relation and exploring whether or not simple 
adjustments in the adopted SNe yields will allow the simu- 
lations to better reproduce the gas-phase silicon abundance 
and stellar metallicities. We leave this for future work. 

Recently, Dave et al. (2008) argued that it is possible 
to reproduce the entropy and metallicity distributions of the 
hot gas in groups with simulations that invoke galactic winds 
with scalings similar to those expected for outflows driven 
by radiation pressure from massive stars on dust grains, but 
which do not include AGN feedback. While it is not straight- 
forward to directly compare the OWLS simulations with 
those of Dave et al. (as there are differences in the imple- 
mentations of star formation, chemodynamics, cooling, and 
feedback), some of their results do not appear to be incon- 
sistent with ours, as we have also demonstrated that models 
without AGN feedback can reproduce these quantities. How- 
ever, their model also appears to yield stellar mass fractions 
in reasonable agreement with the observations, whereas the 
OWLS REF run does not. Plausible reasons for this differ- 
ence are that the energy injected in the momentum-driven 
wind implementation of Dave et al. exceeds that available 
from supernovae (see Schaye et al. 2009). Thus, one must in- 



voke alternative energy sources, namely high energy photons 
exerting pressure on dust grains, to motivate this model. By 
contrast, supermassive BHs are directly observed to be heat- 
ing the hot gas in groups and clusters and the efficiency we 
have adopted has been constrained to yield a good match 
to the present-day black hole mass — halo mass relation and 
the cosmic black hole density. In addition, as Dave et al. 
point out, momentum-driven winds are incapable of shut- 
ting down cooling flows in the most massive haloes at the 
present-day, and therefore yield CBG star formation rates 
well in excess of what is observed. Furthermore, this model 
does not appear to reproduce the large scatter seen in the 
X-ray luminosity-temperature relation of groups. For these 
reasons, we argue that feedback from supermassive BHs is 
likely a key ingredient in galaxy group formation and evo- 
lution. This certainly does not preclude, however, feedback 
from BHs and some form of momentum-driven winds work- 
ing in tandem to shape the properties of groups. 

Additional insights from observations of galaxy groups 
would be very helpful for discriminating between contend- 
ing models and to further constrain the parameters of the 
subgrid physics. In particular, extending entropy measure- 
ments down to lower system masses, where the differences 
between the models should be largest (e.g., Fig. 2), would be 
very useful. Going beyond globally-averaged properties and 
azimuthally-averaged profiles and making detailed compar- 
isons with the full 2D maps of groups (e.g., Finoguenov et 
al. 2006) also represents a very promising avenue for distin- 
guishing between AGN heating mechanisms (e.g., by looking 
at the structure of cavities and sound-wave ripples). 

In addition, adopting well-defined group selection func- 
tions for observed samples that can be applied directly to 
simulated catalogs would be helpful. At present, most sam- 
ples of groups with high quality X-ray observations do not 
have a well-defined selection function. They are comprised 
of systems that have high enough surface brightnesses to 
extract radial profiles, but it is presently not clear whether 
these are 'typical' systems. In the case of massive galaxy 
clusters, for example, the central entropy varies by more 
than an order of magnitude at fixed halo mass (Cavagnolo 
et al. 2009). The systems with the highest central entropies 
(longest cooling times) are referred to as non-cool core sys- 
tems (or non-cooling flow systems, in the old jargon), and 
may represent more than 50% of the cluster population. The 
analogues of non-cool core clusters on the group scale have 
not yet been discovered, maybe because they do not exist, 
or perhaps, more interestingly, because their X-ray surface 
brightnesses are so low that they have escaped detection. If 
a population of X-ray dark groups' could be convincingly 
demonstrated to exist (e.g., in optically-selected group sam- 
ples with follow-up X-ray observations), this would present 
a very interesting challenge to current cosmological simula- 
tions. 
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APPENDIX: EFFECTS OF NUMERICAL 
RESOLUTION 

Here we explore the sensitivity of our results to numerical 
resolution. The full production runs analysed in the main 
body of the paper contain 2 x 512 3 particles, yielding an (ini- 
tial) baryon particle mass of 8.65 x 10 7 ft _1 M Q and a dark 
matter particle mass of 4.06 x 10 s /i" 1 Mq (note that due to 
stellar evolution the mass of gas/star particles can change 
during the course of a simulation, but typically only by a 
small amount). We have also run simulations at lower reso- 
lution using 2 x 256 3 particles, implying the particle masses 
are increased by a factor of 8 over that in the full production 
run. 

In Fig. 13 we compare the median entropy profiles for 
clusters with masses in the range 13.25 ^ log 10 (Msoo /Mq) ^ 
14.25 and the K-band luminosity - uncorrected mean tem- 
perature relations for the full production and lower resolu- 
tion REF runs. The results are robust to a decrease in the 
mass resolution (and we note that the small differences in the 
entropy profiles is well within the system-to-system scatter 
in the entropy profiles). This indicates that the resolution 
of the production runs is adequate to robustly predict the 
properties of galaxy groups. 

This paper has been typeset from a TjtX/ DTffjX file prepared 
by the author. 
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Figure 13. The effect of numerical resolution on the median entropy profile and K-band luminosity - uncorrected mean temperature 
relation for the REF run. (See Fig. 1 and Fig. 7 for comparison.) Overall, the agreement between the two runs is excellent. 



